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Abstract: We present a new technique that combines off-axis Digital 
Holography and Dark Field Microscopy to track lOOnm gold particles 
diffusing in water We show that a single hologram is sufficient to localize 
several particles in a thick sample with a localization accuracy independent 
of the particle position. From our measurements we reconstruct the trajec- 
tories of the particles and derive their 3D diffusion coefficient. Our results 
pave the way for quantitative studies of the motion of single nanoparticle in 
complex media. 
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1. Introduction 

The study of cellular processes at the single-molecule level is a flourishing field of research in 
Biology. Individual molecules labeled with sub-micron markers can now be tracked in a cellu- 
lar environment, and quantitative information about their dynamics can be obtained by recon- 
structing their trajectory. One of the most used techniques for this purpose is single-molecule 
fluorescence microscopy (SMFM), which relies on a labeling with nanometer-sized fluorescent 
markers such as organic dyes or quantum dots. But standard SMFM provides no information on 
the axial position of the marker, limiting this technique to 2D tracking. Recent improvements of 
SMFM such as astigmatism optic [ 1 1 , 4Pi microscopy 1 2 J , double-helix PSF ^ , or multi-planes 
detection BUS) have made possible 3D tracking. Since the depth of field of these techniques is 



limited to a few microns, 3D tracking of molecules that explore larger distances in the thickness 
of a sample requires to continuously adjust the position of the focal plane of the microscope ob- 
jective, which strongly limits the time resolution. Digital Holographic Microscopy (DHM) II6I7I 
circumvents this drawback. In DHM, a CCD camera records the interference pattern between 
the light scattered by the sample and a reference wave, and a single shot is sufficient to deter- 
mine the 3D positions of scatterers embedded in a non-diffusing environment, over a depth of 
typically a hundred of microns. 

As the scattering cross section of a particle scales as the sixth power of its radius HI, how 
easily and accurately a particle can be detected strongly depends on its size. Several publica- 
tions demonstrate the tracking of micron-sized colloids by using in-line holography ll9 tiT5]| . 
with a localization accuracy in the nanometer range through the use of high Numerical Aper- 
ture (NA) microscope objectives. For example, with a lOOx NA 1.4 oil immersion objective, 
Cheong et al. lITOl reported lateral and axial localization accuracies of 4 and 20 nm respectively. 
This result was obtained with polystyrene spheres of diameter d = 1.5 /im, whose scattering 
cross section is quite large. The tracking of li < 100 nm particles, whose scattering cross sec- 
tion is extremely low, is much more difficult, and, as far as we know, has not been demonstrated 
using in-line holography. Yet the detection of such small particles is possible using DHM in a 
off-axis geometry 11611171 with a noise level as low as possible 1181 and using good light scat- 
terers, such as gold nanobeads lfT9l . By this way, = 50 to 200 nm gold particles embedded 
in an agarose gel have been detected and localized ll20l . Since gold particles are not toxic for 
cells, they can be used as markers in biology ||2T1 . and d = 40 nm gold nanobeads fixed on the 
membrane receptors of a living cell have been localized ||22| . More recently, 3D tracking of 
BaTi03 particles with second harmonic generation DHM was also demonstrated 1231. 

Here, the main advantage of DHM, with respect to in-line holography, is the possibility to 
independently adjust the intensity of the illumination and reference beams in order to get the 
best detection sensitivity, by adjusting the reference beam intensity ll24l . and the largest sig- 
nal, by adjusting the illumination beam intensity. Combining DHM with dark field illumination 
allows then to detect nanometer-sized particles, as the sample can be illuminated with an inten- 
sity as large as possible, while avoiding a saturation of the camera 1201 . For example, Atlan et 
al. and Warnasooriya et al. uses Total Internal Reflection (TIR) to detect and localize = 50 
nm and d = 40 nm particles Il20ll22l . But the TIR configuration used in these experiments 
yield a standing wave which does not allow to track moving particles : when a moving particle 
crosses a node, the illumination (and thus the signal) goes down to zero, and the particle is lost. 
Dubois et al. uses another dark field illumination configuration that focuses the illumination 
on a mask f25\. Since the illumination is parallel to the optical axis, no standing wave can ap- 
pear, but since the illumination passes through the microscope objective, one expects parasitic 
reflections of the illumination beam. 

In this paper, we present a Digital Holographic Microscopy technique which makes possible 
to track = 100 nm gold particles 3D diffusing in water The illumination is parallel to the 
optical axis to prevent the formation of standing waves, and the holographic signal is collected 
by a NA=0.5 dark field reflecting microscope objective. With such objective, the illumination 
beam is masked before the microscope objective and parasitic reflections are avoided. This 
yields high dynamic dark field illumination, which makes possible to detect, localize and track 
d = 100 nm particles. First we describe the setup, which combines dark-field microscopy and 
off-axis holography. Then we present the algorithm of reconstruction, our procedure to localize 
the beads, and describe how we can reach a real-time localization by performing calculations on 
a graphic card. Finally we show that our setup allows us to track gold nanoparticles in motion 
with a lateral {x,y) resolution of 3 nm and an axial (z) resolution of ~ 70 nm. Since NA=0.5, 
the resolution (especially in z) is lower than with NA=1.4 in-line holography lITOl . We also 



show that the depth of field of our holographic microscope is made two orders of magnitude 
larger than in optical microscopy. 

2. Digital holography setup 

Our DHM experimental setup, depicted in Fig.[Tl is designed to investigate the Brownian mo- 
tion of lOOnm diameter gold particles diffusing in a 17 x 3.8 x 0.4mm (length x width x 
height) cell chamber (Ibidi ©/X— Slide) filled with water The concentration of nanoparticles is 
adjusted to 5.6 10^ particles/mm^ to have a few particles per field of view. 

The light source is a 660 nm Diode Pump Solid State laser (Crystal Laser ©) with a short co- 
herence length (^ 600 jim) to avoid parasitic interferences raising from reflexions between the 
optical elements. The laser beam is split into two beams by a Polarizing Beam Splitter (PBS), 
a half-wave plate before the PBS setting the ratio of energy between the emerging beams. The 
reference beam passes through a dove prism fixed on a micrometer translation stage to adjust 
the length of the optical path. This beam is spatially filtered through a 35jLim diameter pinhole 
and then expanded as to uniformly cover the CCD chip of the camera (512 x 512 pixels, Andor 
Luca R ©). The illumination beam is focalized on the sample with a plano-convex lens of focal 
length 12.5 cm (waist diameter ^ 200 /xm, laser intensity ^ 250 W/cm^). The light scattered 
by the beads is collected in transmission with a dark-field reflecting objective (Edmund Optics 
©, ReflX series) of NA = 0.5 and 36X magnification. A small mask on the input of the objec- 
tive limits the collection of light between NA = 0.2 and NA = 0.5, so the illumination beam is 
totally blocked after passing through the sample. This dark-field configuration prevents the sat- 
uration of the CCD chip. A non-polarizing 50/50 beam splitter behind the objective combines 
the scattered light with the reference beam. The CCD camera records the interference pattern 
on 16bits frames with a 22.5 Hz rate. The last beam-splitter is tilted by few degrees to be in 
off-axis configuration. 



Laser PBS 50/50 ND filter A/2 wave plate 

@660 nm 

Fig. 1. Experimental setup. The sample is located in the X,Y plane. Z is the optical axis of 
the microscope objective 

3. Numerical reconstruction 

3.1. Reconstruction of the scattering field 

In digital holography, the CCD sensor records an intensity Iced which is the interference of the 
reference beam with the light scattered by the nano-objects. I^cd is thus a sum of 4 terms: 
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(1) 



where 7^/ is the intensity of the reference beam, Iscatt the intensity of the scattered Ught and 
Eref,Escatt are the electric fields for the reference beam and the scattered light respectively. As 
the scattering cross section of a lOOnm diameter gold particle calculated with the Rayleigh-Mie 
scattering model is about 0.015 /im^ at 660 nm [Sl, the integration over the collection solid angle 
of the objective (0.2 < NA < 0.5) gives Iscarr — 3 x 10^^ x lutumination for a single particle. Thus 
Iscatt can be neglected compared to the other terms of Eq.([T]). The field Escmt, the amplitude of 
which is about \/ Iscatt, contains the phase of the scattered light necessary to a 3D localization 
of the particle. In the Fourier space, the tilted beam splitter adds a spatial frequency on the two 
conjugates terms of interference, and thus the different terms of Iced are spatially separated, i.e. 
Iref remains centered on the zero frequency of the Fourier plane while the two cross terms of 
interference are centered on the spatial frequency induced by the off-axis geometry. 

Since 3D reconstruction is a time consuming task even using recent multi-cores proces- 
sors, we developed parallel calculations on a Graphic card (Nvidia Geforce GT470, 448 cores) 
ll26]^9l . Our C++ based algorithm uses the Nvidia CUDA library to decompose the calcula- 
tions on the GPUs of the card. Among the existing methods to reconstruct holograms, the most 
common is the convolution method described by Schnars et al. Q. The main drawback of this 
method is that the pixel size of the reconstructed image depends on the distance of reconstruc- 
tion, so that the image of a thick sample is distorted. Therefore this method is not convenient 
for 3D tracking, as the lateral scale depends on the depth of reconstruction considered. Here 
we chose to reconstruct the holograms with the angular spectrum method II30II31I . which, by 
compensating the sphericity of the signal wave, allows to reconstruct the hologram without 
distortion. Our algorithm can be decomposed in 6 steps : 

i. subtraction of the background: in order to increase the signal-to-noise ratio, we subtract 

from the last recorded frame the average of the ten previous frames. Phase shifting holog- 
raphy 1321 is also an effective technique for reducing noise, but the minimal delay Af = 44 
ms between two frames, which is driven by our camera, is too long to use this technique 
for nanoparticle tracking. This step of calculation is performed only for particles in mo- 
tion (i.e. we skipped this step for the results presented in |4.1| i. 

ii. numerical correction of the signal wave sphericity: the hologram I^cd is multiplied by a 

complex phase matrix M to compensate the sphericity induced by the microscope ob- 
jective on the signal wave: 



where d is the local radius of curvature of the wave on the CCD plane. If the reference 
wave is a plane wave, this distance d is also the distance between the CCD chip and the 
back focal plane of the objective. 

m. first FFT: the direct Fourier transform of the corrected hologram is calculated using the 
CUDA CUFFT library: 



Fig.|2a) shows the intensity |^| in the k-space, in logarithmic scale in the case where the 
background is not removed (step (i) skipped). In the middle of Fig. I2a), the zero-order 
appears as a square because of the multiplication by the matrix M. The term related to 
Escatt is in the down-right corner, centered on the spatial frequency induced by the off- 
axis geometry. The term related to E*^^^ is centered on the conjugate frequency (top-left 
corner). At this step, the calculation is equivalent of the reconstruction in one FFT (con- 
volution method) of the hologram at the distance d described in the previous step. Since 




(2) 




(3) 



the back focal plane of our microscope objective coincides with the output pupil plane, 
which is common for high magnification objective, we see on the down-right corner a 
sharp reconstruction of the output pupil plane. If we change in step (i) the parameter d to 
—d, the image of the output pupil would be sharp in the top-left corner, while the term 
related to Escatt would be blurred. 

Fig-Etb) shows \H\ in the k-space when the background is removed (step (i) performed). 
The zero-order term in the middle of the k-space is largely removed compare to Fig.|2a), 
which reduces the recovery between the zero-order and Escatt- 

iv. spatial filtering and centering: to remove the zero-order term and replace the term related 
to Escatt in the middle of the Fourier plane, a round numerical filter which matches with 
the output pupil of the objective is applied. Since the shape of the pupil is sharp in the 
k-space, we can isolate precisely the pixels containing the signal, minimizing the loss of 
information. To more precisely calibrate the radius and the center of the filter, we used 
a diffusive paper as a sample before performing experiments. Since the paper scatters 
light uniformly, all the spatial frequencies that the microscope objective can collect are 
recorded. Fig. |2jd) shows the intensity |^| in logarithmic scale when the sample is re- 
placed by a diffusing paper We clearly see the shape of the output pupil of the objective 
and the dark-field mask in the center. We set the filter mask to match with the shape of 
the output pupil (white dotted circle in Fig.|2d)). The filtered part is then translated into 
the middle of a 512 x 512 calculation grid in order to compensate the off-axis shift. 

y. propagation: the matrix obtained is multiplied by a propagation matrix K{kx,ky,z), which 
is the exact form of the matrix propagation as given by Kim et al. [i33l : 

K{k^,ky,z) = exp (iz X ^kl-kj-kj^ (4) 

where 

2mi 2n{x-256) 2n{y-256) 
h = ^-, k^ = -— — , ky^—— — (5) 

to propagate the hologram by a distance z in the axial direction. Ap,-, is the magnified 
pixel size. These equations are suited for holograms of 512 x 512 pixels. 

vi. second FFT: finally the inverse FFT is calculated. 

For each hologram, the steps (v) and (vi) are repeated in order to get a stack of the scattered 
field at different depths, with a propagation step 5z = lOOnm: 

H{x,y,n-8z) = FFT ' [H{k„ky) x K{k,,ky,n-dz)] , (6) 

where n is an integer 

An example of reconstruction is shown in Fig.|2|c). We can see a nanoparticle (1) localized in 
the considered plane and the intensity of the field scattered by two particles at different depths (2 
and 3). The reconstruction of a 50 jUm thick volume requires to calculate one FFT to reconstruct 
the hologram in the output pupil plane, then to calculate 500 inverse FFTs for the slices of the 
stack. These 501 FFTs typically require one minute when the calculation is performed on the 
CPU, even with a recent multi-cores computer We reduce this time by a factor of 30 when the 
calculation is parallelized on the graphic card. 




Fig. 2. a: Intensity \H{kx,ky)\ in logarithmic scale. The zero-order term appears as a red 
square distorted by the multiplication with the phase matrix M. The term of interest is 
located in the down-right corner b: \H(kx,ky)\ when the average of the ten previous holo- 
grams is subtracted before calculating the FFT The zero-order term is largely removed, so 
the recovery between this term and the region of interest in the down-right comer is re- 
duced, c: Two-dimensional reconstruction at a fixed depth of the sample. A gold nanobead 
(1) is localized in this plan and the shape of the intensity of the field scattered by other 
beads at other depths (2 and 3) is visible, d: Intensity \H{kji,ky)\ in logarithmic scale when 
the sample is replaced by a diffusive paper. The area of the output pupil of the objective is 
sharply defined. The white-dotted circle shows the mask of the numerical filter used for the 
reconstruction. 

3.2. Method of localization 

Once the three dimensional map of the scattered field is calculated, the beads are localized by 
pointing the local maxima of the field's intensity. Fig. [3] shows the field's intensity in the X,Y 
and Z directions. The Full Width at Half Maximum (FWHM) is about 1 jttm in lateral direction, 
and 14/zm in the axial direction as expected for the Point Spread Function associated with 
our microscope objective II34II35II . In order to reach a sub-pixel size resolution, the intensity 
is fitted by a Gaussian curve in X and Y using the pixel for which the intensity is maximum 
and the two adjacent pixels. As the intensity profile in Z is not gaussian and 10 times larger 
than in X and Y, we fit the maximum of the peak with a parabola using the pixel / for which 
the intensity is maximum and the two adjacent pixels ; — 2 and ; + 2 (Fig. Od)). We chose this 
simple localization method because programming a more elaborate fit, as T-Matrix theory based 
computation 151, with CUDA (programming language of the graphic card) is more complicated 
and would considerably slow down the process. Our method shows good performance (see Fig. 



|4]and Fig.O, but a better resolution may be achieved using T-Matrix theory. 
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Fig. 3. PSF of the field's intensity for a lOOnm nanobead, (a) in X, (b) Y, (c) and in the 
Z direction, (d) shows the top of the curve (c) (blue points) and the polynomial fit (green 
solid line). 



4. Results 

4.1. Embedded particles in a gel 

To evaluate the lateral and axial localization accuracy of our setup, we localize a single lOOnm 
diameter gold particle embedded in a 1% agarose gel. The localization accuracy is evaluated 
by calculating the standard deviation of 200 positions of the bead obtained from successive 
frames, with an exposure time for each frame 5f = 1 ms. For a particle in the focal plane of the 
objective, we found a lateral localization accuracy of 3 nm and an axial localization accuracy 
of ~ 150nm (Fig. EJa) andUb)). Then the distance between the particle and the focal plane 
was increased by steps of 20 /im. For each step we recorded 200 holograms and determined the 
mean position of the bead as well as the lateral and axial accuracy. The Fig. |5] compares the 
mean axial position of the particle with the mechanical displacement along Z of the sample, 
which are in excellent agreement. The localization accuracy in X,Y and Z as a function of the 



axial position of the bead is shown in Fig.|Ha) andHfb). While the lateral accuracy is constant 
(~ 3nm) for |Z| < 250 nm, the axial accuracy slightly depends on Z. It is about ~ 150nm 
around Z = 0, then decreases to ^ 70 nm for |Z| < 250 nm. This accuracy strongly increases 
for |Z| > 250/xm, and for Z > 400 /xm, the localization of the particle is not possible because 
the scattered signal level reaches the noise level. The local maximum at Z = observed on the 
axial accuracy curve (Fig.Ub)) shows that the localization is not optimal when the gold particle 
is the focal plane of the objective. In this case the particle is imaged on a small area of the CCD 
chip, so that the interference pattern spreads on a small number of pixels, which degrades the 
quality of the reconstruction ll36l . 

In the case of particle tracking, the position of the focal plane in the sample has to be fixed. To 
minimize the spherical aberrations due to the presence of the coverslip, the focal plane should 
coincide with the sample-coverslip interface. In this configuration, a moving particle cannot 
cross the focal plane, so that the localization accuracy remains optimal and constant over a 
depth of 250 /im, that corresponds to the part between abscissa and abscissa 250 in Fig. ^a.) 
andUb). Yet tracking particles remains possible when the focal plane is set above the sample- 
coverslip interface as demonstrated below, as the cost of a slightly worse axial locaHzation 
accuracy. 




4.2. Particles in Brownian motion 

We now consider gold particles in Brownian motion. According to the Stokes-Einstein equa- 
tion, the diffusion coefficient D of the particles is given by: 

D= =4.2±0.2^m-s"', (7) 
onrjr 

where is the Boltzmann constant, T the room temperature (20 °C), rj the viscosity of water 
(1.0 mPa.s at 20 °C) and r = 50 ± 2 nm the radius of the nanobead (size dispersity given by the 
provider BBInternational). The exposure time is 5f = 1 ms and the time between two frames 
At = 44 ms. The mean distance traveled along one direction by a Brownian particle during 8t is 
90 nm, which is smaller than the lateral size of the magnified pixels (160 nm). Consequently, the 
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Fig. 5. Mean axial position of an embedded particle, estimated using the reconstruction 
algorithm, as a function of the mechanical displacement of the sample along Z. 

signal from a particle is not blurred over several pixels during St. The mean distance covered 
along a given axis during At is 620nm, which corresponds to approximately 5 pixels. 

As shown in Fig.|6] our method allows us to simultaneously track several particles. A volume 
of 80 X 80 X 250jUm(X xY xZ), i.e. 512 x 512 x 400 pixels, can be reconstructed from a single 
hologram, and the localization method described in the section fOl can be performed to localize 
several beads with a sub-pixel accuracy. By repeating the algorithm for successive frames, we 
could for instance reconstruct the trajectories of 3 gold particles diffusing in water (Fig. |6l). 
We were able to track particles during up to 10 s 200 frames). Since the time needed to 
reconstruct a volume of 512 x 512 x 400 pixels is about 0.5 s (i.e. much larger than At ~ 44 
ms), reconstruction is necessarily a post-processing process in this case. 

In contrast, when a single particle is tracked, 3D localization in a given frame requires to 
reconstruct only a few slices around the position of the particle in the previous frame. For a 
lOOnm particle, the mean distance traveled along Z during At is ~ 600 nm, thus only 24 recon- 
structions (12 reconstructions above/below the previous position of the particle) are sufficient 
to find the new position of the particle. The calculation of these 24 slices from the hologram re- 
quires 33 ms, which is smaller than At = 44 ms. Real-time tracking is thus possible for a single 
particle under the condition that reconstruction is performed fast enough, for instance by using 
a graphic card unit as described above. 

To evidence the Brownian motion, we calculated the 3D Mean Square Displacement (MSD) 
of a bead from the red trajectory plotted in Fig. |6l 



where = 200 is the total number of positions, as well as the ID MSD along the directions X, 
Y, and Z. For a Brownian motion with a diffusion constant D, the MSD curve depends linearly 
on time, and the slope of the curve is 2"'' D, where n^/ = 3 for the 3D MSD (Fig. |2ta)) and 
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Fig. 6. 3D trajectories of 3 particles (red, green, blue) reconstructed from 200 successive 
frames. The focal plane (Z = 0) was set at about 150 flm above the coverslip. Although this 
setting is not optimal for tracking a particle diffusing around the focal plane, as explained 
in section I4J] trajectory reconstruction is still possible (red trajectory). 

«(/ = 1 for ID MSD (Fig.|7tb), (c) and (d)). As expected for a Brownian motion, experimental 
3D and ID MSD depend linearly on time, and a fit of the first six points of the curves gives 
D^4.3±0.5lJ.m^s-^ along X, D = 5 .0 ± 0.3 ^m^s" ' along F, D = 2.9 ± 0.5 ^m^s" ' along Z 
and D = 4. 1 ± 0.5 jUm^s" ^ for the 3D motion. These values are in agreement with the theoretical 
value D = 4.2 ± 0.2 /im^s" ' predicted by Eq. 

5. Conclusion 

In this paper, we show that Digital Holographic Microscopy can be used to track d = lOOnm 
gold particles diffusing in water As the intensity of the light scattered by such nanoparticles is 
height orders of magnitude smaller than that of the excitation light, we combined holography to 
standing wave free dark-field microscopy to completely block the illumination beam, thereby 
preventing a saturation of the CCD chip of the camera. A single hologram, recorded with an 
exposure time of only 1 ms, is sufficient to localize several particles in a 250 /im thick sample, 
with a lateral {x,y) localization accuracy of ~ 3nm and an axial (z) localization accuracy of 
~ 70 nm. As our dark field microscope involves a NA=0.5 reflecting microscope objective, the 
resolution, especially in z, obtained here with nanoparticles (c/ = 100 nm), is lower than the 
resolution reached in in-line holography using a NA=1.4 objective and micron-sized objects 
(d ~ 1.5 /im) flOll. This is the cost to pay for the detection of nanometer-sized particles. We 
were able to reconstruct particle trajectories, evidence the Brownian nature of the motion and 
determine the related diffusion coefficient. The accuracy achieved by our setup is comparable 
with that reached with super-resolution microscopy: a localization accuracy of 20 nm in X 
and Y, and 50 nm in Z has been reported for fluorescence dye or quantum dots imaged with 
STORM microscopy using optical astigmatism HI, and an accuracy of 10 nm over a depth 
of 2 jim has been reached using a double helix PSF |f3l| . Severals results of 3D tracking of 
quantum dots using PALM microscopy report an axial accuracy between 10 and 75 nm, over 
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Fig. 7. Mean Square Displacement curves derived from the trajectory of a nanoparticle in 
Brownian motion (200 points, time step A/ = 44 ms). (a) 3D MSD, (b) MSD along X, (c) 
MSD along Y and (d) MSD along Z. Blue line: linear fit over the first 6 points of the MSD. 

a depth of typically 1 jxm ||4]|5]. PALM and STORM microscopy reach higher 3D localization 
accuracies than our DHM setup, but over a depth two order of magnitude lower. This ability to 
track nanoparticles up to 250/im from the focal plane with a constant localization accuracy is 
the strength of DHM compare to super-resolution techniques. Our results pave the way for the 
use of gold nanobeads as markers in more complex media such as cellular environment. 
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